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Abstract 

We report on a methodology for the evaluation of the DC characteristics, small-signal frequency 
response and large-signal dynamic response of carrier and photon density responses in semicon- 
ductor laser diodes. A single mode laser is considered and described with a pair of rate equations 
containing a novel non-linear gain compensation term depending on a single parameter that can 
be chosen arbitrarily. This approach can be applied to any type of solid-state laser as long as it is 
described by a set of rate equations. 
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I. INTRODUCTION 



Generally, lasers are described with a set of rate equations describing the generation- 
recombination of carriers or emission-absorption of photons. The equations are typically 
systems of first-order differential equations belonging to the population evolution type. 

In this report we specialize to the simple case of single-mode laser diode where the 
system of equations reduces to two: one for the carrier population (density) and another for 
the photon population (density). We concentrate on the case the equations embody a novel 
non-linear gain compensation term that depends on a single parameter that can be chosen 
arbitrarily. The effect of this parameter is studied on the DC characteristics, small-signal 
frequency response and large-signal dynamic response of the laser. 

The Single Mode (SM) laser diode rate equations are written as: 



N represents the electron density (N t at transparency) and 5* the photon density. r n 
is the electron spontaneous lifetime and r p h is the photon lifetime. (3 is the fraction of 
spontaneous emission coupled to the lasing mode, G the optical confinement factor, g is the 
differential gain and e is the gain compression parameter, q is the electron charge, V the vol- 
ume of the active region and / is the injection current which, in general, is a function of time. 

The main novelty in these equations lie in the non-linear gain factor that has been 



traditionally modelled by the term g(l — eS), g/(l + eS) or even g/J(l + eS) |TJ. In this 
study, this term is taken to an arbitrary power b in the interval [0,1.5]. In order to cover 
in the simulation the case g(l — eS), we formally allocate b=-l to this case. We analyse 
in this work the effects b has on the static, small- signal frequency response as well as the 
large-signal temporal variation of the carrier and photon densities. 

This report is organized as follows: in section 2 we outline the evaluation of the Laser 
S — I (photon density versus injection current) DC characteristics. In section 3, the small- 
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signal frequency response is derived and in section 4, we illustrate the laser response to 
time-dependent injection currents and section 5 contains our conclusions. 



II. DC CHARACTERISTICS 



The DC limit of the SM laser rate equations is given by: 



/ N 

-T7---9(N-N t ) 

qV r n 



(1 + eS) c 



= 



T/3N 



— + r g (N 

Tph 



N 



S 



1 + eS) 



(3) 
(4) 



From (4) we extract the value of (N — N t ) as: 
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and substitute it in (3) relating I to S directly. This results in: 
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This relation gives / as a function of S. Generating a series of data values (S, I) by 
varying S and reverse writing them as (J, S) will result in the S — I DC characteristic. 



We consider two laser models A and B (Physical parameters given in Appendix A) and 
illustrate this procedure with the characteristics displayed in Figures 1 and 2 for different 
values of b. It appears from the figures that b does not affect significantly the DC charac- 
teristics in contrast with the frequency response and dynamic response as seen in the next 
sections. 



III. SMALL-SIGNAL FREQUENCY RESPONSE 

In order to derive the small signal frequency response, we assume all quantities /, iV and 
S are taken around some equilibrium values Jq, N q and So and hence: 



I = I + 6I(t), N = N + 5N(t), S = S + SS(t) 
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(7) 



This means that equation (1) under variation (7) reads: 
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In order to tackle the small-signal frequency response we switch to the time-harmonic case 
where the time derivatives are given by: dSI/dt = juj5I,d5N/dt = ju5N,d5S/dt = juSS. 
This results in a system of equations relating the three variations 51, 5N and 5S: 
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Taking the ratio of the above yields the small-signal frequency response: 
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where E, G and C are given by: 
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The standard normalised form ( dB at frequency) of the frequency response is taken 



as: 
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It is displayed for laser A and B in Figures 3 and 4 for different values of b. 

In sharp contrast with the DC case, the value of b, has a drastic effect on the frequency 
response as observed in the above figures. We differ from Way |3| in the frequency response 
due to a discrepancy in the estimation of the resonance frequency. Way defines the resonance 
frequency as (using our own notation and adapting it to Way's [§J case): 



where / is the bias and Ith is the threshold current. Our corresponding formula by 
inspection of (16) would be: 



So that are found numerically with a Newton method |TT[ adapted to (5) and (6). 



(17) and (18) since he wanted to get an analytical estimate of the resonance frequency. In 
order to estimate the discrepancies between our work and Way's we display in Figure 5 the 
frequency response for various values of the bias current expressed in terms of the threshold 
current Ith ~21mA. 

IV. LARGE-SIGNAL DYNAMIC RESPONSE 

We exploit three possible excitation scheme for the injection current: 

1. A step excitation in order to evaluate the step response of the laser. 

2. A graded response with a Gaussian time excitation towards a higher injection level. 

3. A modulation injection in order to estimate the modulation response of the laser for 
large excursions of the injection current. 




(17) 




Obviously, Way's estimation [[| of the resonance frequency is approximate by comparing 
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As an example we consider lasers A and B biased at t=0 and excited with an additional 
square pulse triggered after 5 nanoseconds of operation. Figures 6-9 show the large signal 
dynamic responses for the carrier and photon densities as functions of time for different 
values of b. As expected, the value of b deeply affects the dynamic response of the Laser as 
observed in the figures 6-9. 



V. CONCLUSIONS 



We have developed an approach that evaluates the DC characteristics, small-signal 
frequency response and large- signal dynamic responses of carrier and photon densities 
in single-mode semiconductor laser diodes directly from the rate equations. The laser is 
described with a pair of rate equations containing a novel non-linear gain compensation 
term depending on a single parameter b that can be chosen arbitrarily in the range [0,1.5] 
as has been shown recently with a microscopic calculation of plasma-heating induced 
intensity-dependent gain effects [Q]. 

Our DC evaluations agree with several published results [@]> ^ U, U H and our 
large-signal dynamical reponses agree also with what is well established in the literature 
0. We differ with the small-signal frequency response results given by Way || due to a 
discrepancy in the frequency response and resonance frequency estimation. 



The parameter b has almost no effect on the DC characteristics but deeply affects the 
small-signal frequency response as well as the large-signal dynamic response. 

The methodology developed in this work can be easily generalised to an arbitrary number 
of state equations that appear in multimode semiconductor lasers, MQW (Muli-Quantum- 
Well) or Strained layer lasers |p], p], pTO . 
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Appendix 



LASER A (ORTEL LS-620) 


G (mode confinement factor) 


0.646 


r n (electron spontaneous lifetime) 


3.72 10~ 9 sec 


T p h (photon lifetime) 


2.10- 12 sec 


N t (electron density at transparency) 


4.6 10 24 m- 3 


g (differential gain ) 


10~ 12 m 3 /sec 


e (gain compression parameter) 


3.8 10~ 23 m 3 


P (fraction of spontaneous emission coupled to the lasing mode) 


0.001 


V (volume of the active region) 


9.10~ 17 m 3 


TABLE I: Physical parameters of Laser A 


LASER B 


G (mode confinement factor) 


0.34 


T n (electron spontaneous lifetime) 


3. 10~ 9 sec 


Tph (photon lifetime) 


2.10- 12 sec 


N t (electron density at transparency) 


10 24 m- 3 


g (differential gain ) 


3.10~ 12 m 3 /sec 


e (gain compression parameter) 


3. 10~ 23 m 3 


(5 (fraction of spontaneous emission coupled to the lasing mode) 


0.001 


V (volume of the active region) 


3.6 10" 18 m 3 



TABLE II: Physical parameters of Laser B 
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Figure Captions 



Fig. 1: S — I characteristics (Laser A) for various values of b: -1.0, 0.5, 1.0 and 1.5. The 
different values of b do not affect significantly the DC characteristics. 

Fig. 2: S — I characteristics (Laser B) for various values of b: -1.0, 0.5, 1.0 and 1.5. The 
different values of b do not affect significantly the DC characteristics. 

Fig. 3: Small signal frequency response amplitude (Laser A) versus frequency for various 
values of b: -1.0, 0.5, 1.0 and 1.5. The bias current is chosen to be 40 mA. 

Fig. 4: Small signal frequency response amplitude (Laser B) versus frequency for various val- 
ues of b: -1.0, 0.5, 1.0 and 1.5. The bias current is chosen to be 1 mA. 

Fig. 5: Small signal frequency response amplitude (Laser A or Way's case [4]) versus frequency 
for various values of the bias current. The bias current is taken as I t h, 1.25 I t h, 1-75 
Ith and 2.5 I t h where I t h ~ 21mA. 

Fig. 6: Laser A large signal dynamic response amplitude of the carrier density versus time 
for various values ofb: -1.0, 0.5, 1.0. The bias current is 40 mA and a square pulse 
excitation of 10 mA is applied after 5 nanoseconds. The dynamic response for b=1.5 
is off the graph scale. 

Fig. 7: Laser A large signal dynamic response amplitude of the photon density versus time 
for various values ofb: -1.0, 0.5, 1.0. The bias current is 40 mA and a square pulse 
excitation of 10 mA is applied after 5 nanoseconds. The dynamic response for b=1.5 
is off the graph scale. 

Fig. 8: Laser B large signal dynamic response amplitude of the carrier density versus time 
for various values ofb: -1.0, 0.5, 1.0. The bias current is 1 mA and a square pulse 
excitation of 0.5 mA is applied after 5 nanoseconds. The dynamic response for b=1.5 
is off the graph scale. 

Fig. 9: Laser B large signal dynamic response amplitude of the photon density versus time 
for various values ofb: -1.0, 0.5, 1.0. The bias current is 1 mA and a square pulse 
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excitation of 0.5 mA is applied after 5 nanoseconds. The dynamic response for b=1.5 
is off the graph scale. 
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FIG. 1: S — I characteristics (Laser A) for various values of b: -1.0, 0.5, 1.0 and 1.5. The different 
values of b do not affect significantly the DC characteristics. 
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FIG. 2: S — I characteristics (Laser B) for various values of b: -1.0, 0.5, 1.0 and 1.5. The different 
values of b do not affect significantly the DC characteristics. 
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FIG. 3: Small signal frequency response amplitude (Laser A) versus frequency for various values 
of b. From top to bottom, b is -1.0, 0.5, 1.0 and 1.5 respectively. The bias current is chosen as 40 
mA. 
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FIG. 4: Small signal frequency response amplitude (Laser B) versus frequency for various values 
of b. From top to bottom, b is -1.0, 0.5, 1.0 and 1.5 respectively. The bias current is chosen as 1 
mA. 
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FIG. 5: Small signal frequency response amplitude (Laser A or Way's case [4]) versus frequency 
for various values of the bias current. From top to bottom, the bias current is taken as I t h, 1-25 
Ith, 1-75 Ith and 2.5 I th where I th ~ 21mA. 
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FIG. 6: Laser A large signal dynamic response amplitude of the carrier density versus time for 
various values of b: -1.0, 0.5, 1.0. The bias current is 40 mA and a square pulse excitation of 10 
mA is applied after 5 nanoseconds. The dynamic response for b=1.5 is off the graph scale. 
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FIG. 7: Laser A large signal dynamic response amplitude of the photon density versus time for 
various values of b: -1.0, 0.5, 1.0. The bias current is 40 mA and a square pulse excitation of 10 
mA is applied after 5 nanoseconds. The dynamic response for b=1.5 is off the graph scale. 
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FIG. 8: Laser B large signal dynamic response amplitude of the carrier density versus time for 
various values of b: -1.0, 0.5, 1.0. The bias current is 1 mA and a square pulse excitation of 0.5 
mA is applied after 5 nanoseconds. The dynamic response for b=1.5 is off the graph scale. 
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FIG. 9: Laser B large signal dynamic response amplitude of the photon density versus time for 
various values of b: -1.0, 0.5, 1.0. The bias current is 1 mA and a square pulse excitation of 0.5 
mA is applied after 5 nanoseconds. The dynamic response for b=1.5 is off the graph scale. 
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